Nonlinear electrokinetics at large applied voltages 
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The classical theory of electrokinetic phenomena assumes a dilute solution of point-like ions in 
chemical equilibrium with a surface whose double-layer voltage is of order the thermal voltage, 
ksT/e = 25 mV. In nonlinear "induced-charge" electrokinetic phenomena, such as AC electro- 
osmosis, several Volts ~ lOOfcsT/e are applied to the double layer, so the theory breaks down 
and cannot explain many observed features. We argue that, under such a large voltage, counterions 
"condense" near the surface, even for dilute bulk solutions. Based on simple models, we predict that 
the double-layer capacitance decreases and the electro-osmotic mobility saturates at large voltages, 
due to steric repulsion and increased viscosity of the condensed layer, respectively. The former 
suffices to explain observed high frequency flow reversal in AC electro-osmosis; the latter leads to 
a salt concentration dependence of induced-charge flows comparable to experiments, although a 
complete theory is still lacking. 
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Electrically driven flows in ionic solutions are finding 
many new applications in microfluidics [l| . The theory of 
electro-osmosis Q was developed for slip past a surface in 
chemical equilibrium, whose double-layer voltage is typ- 
ically of order fc^T/e = 25 mV. However the discover 
of AC electro-osmosis (ACEO) at micro-electrodes 0, 
has shifted attention to a new regime, where the induced 
double-layer voltage is w lOOfcsT/e, oscillating at fre- 
quencies up to 100 kHz, and nonuniform at the micron 
scale. Related effects of induced-charge electro-osmosis 
(iceo) mm also occur around colloidal particles |7fl and 

polarizable microstructures § (in AC or DC fields), and 
asymmetric particles can move by induced-charge elec- 
trophoresis (ICEP) di]. 

In all of these situations, low-vo ltag e theories fail to 
predict crucial experimental trends llfj , such as flow de- 
cay at high salt concentration 11 , [Til, EH and flow re- 
versal in asymmetric pumps at high voltage and high 
frequency [ll, 14, 15j]. In this Letter, we attribute these 
failures to the breakdown of the dilute solution approx- 
imation at large applied voltages. Based on very simple 
models, we predict two general effects due to counterion 
crowding - decay of the double-layer capacitance and sat- 
uration of the electro-osmotic mobility - which begin to 
explain the experimental data. 

Experimental puzzles - ICEO flows are rather complex, 
so many simplifications have been made to arrive at an 
operational model [!, 0, @, 0, 0- For thin double layers, 
the first step is to solve Laplace's equation (Ohm's Law) 
for the electrostatic potential in the conducting bulk, 
V 2 = 0, with a capacitance-like boundary condition 
to close the "RC" circuit Hi, 
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ductivity a and diffuse-layer capacitance Cd are usu- 
ally constants, although these assumptions can be re- 
laxed ■ The second step is to solve for a Stokes flow 
with the Helmholtz-Smoluchowski effective slip boundary 
condition, 



u s = -6E i = -^E t 
Vb 



(2) 



where the local diffuse-layer voltage 4 r rj(0) responds to 
the normal electric field E n — — h ■ V</>; the bulk con- 



where b(^D) is the electro-osmotic mobility, E* the tan- 
gential field, C = is the "zeta potential" at the shear 
plane (£ = ^f> d in the simplest models) , and £f> and rjh are 
the permittivity and viscosity of the bulk solvent. While 
this model based on dilute solution theory can only be 
justified for 'F r> ~ ksT/e [6|, ll6|, Ll8| , it describes many 
features of ICEO flows at moderate and large voltages. 
Nevertheless, some crucial effects are still missing. 

ICEO flows have a strong sensitivity to solution chem- 
istry, which is under-reported and unexplained. Recent 
experiments suggest a universal logarithmic decay of the 
mobility, b oc log(c c /co), with bulk concentration cq seen 
in KC1 for ACEO micropumps [llj, in KC1 and CaCo3 
for ICEO flows around metal posts and in NaCl 

for ICEP of colloids with partial metal coatings In 
all cases, the flow practically vanishes for Co > c c « 30 
mM, well below the concentration of most biological flu- 
ids (cq w 0.3 M). ICEO flows and AC-field induced 
interactions in colloids [l9| are also sensitive to the par- 
ticular ions, at a given concentration. 

The low-voltage model also fails to describe the re- 
versal of ICEO flow observed in some (but not all) ex- 
perimental situations. Flow reversal was first reported 
around metal particles in water [20| . where the velocity 
agreed with the theory Q only for micron-sized particles 
and reversed for larger ones (90 — 400/im). The transition 
occured when several volts was applied across the parti- 
cle and reversal was attributed to Faradaic reactions [13] • 
Flow reversal has also been observed at high voltage (> 2 
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V) and high freq uency (10-100 kHz) in ACEO pumping 
of dilute KC1 [lllllJ] and deionized water 15[ with 10/im 
scale electrode arrays. This reversal was first attributed 
to Faradaic reactions lj], but simulations with Butlcr- 
Volmer reaction kinetics have failed to predict the ob- 
served flow • With non-planar 3D electrodes 21 1 , the 
low-voltage model also fails to predict the double-peaked 
frequency spectrum, which accompanies flow reversal in 
some geometries [l 51 ] . 

Although Faradaic reactions surely occur at large volt- 
ages, they are dominant at low frequencies in ACEO 
simulations and experiments (when gas bubbles 

arise) [ll| . Dilute solution theories also predict that non- 
linear effects dominate at low frequencies: The differen- 
tial capacitance of the diffuse layer Q , 



C D {^! D ) = -^cosh 
Ad 
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causes the RC charging time to grow exponentially with 
voltage [l(J, and salt adsorption and tangential conduc- 
tion by the diffuse layer are coupled to (much slower) 
bulk diffusion [lg, |l7j. Strong flow reversal is experi- 
mentally observed at much higher frequencies, and the 
unexplained concentration dependence seems to be inde- 
pendent of frequency. We conclude that dilute solution 
theories do not properly describe the dynamics of elec- 
trolytes at large voltages. 

Crowding effects - All dilute solution theories, which 
describe point-like ions in a mean-field approximation, 
break down when the crowding of ions become significant, 
when steric effects and correlations potentially become 
important. If this can be translated into a characteristic 
length scale a for the distance between ions, then the 
validity of the Poisson-Boltzmann (PB) is limited by a 
cutoff concentration c max — a~ 3 , reached at a fairly small 
voltage, 



H" 3 c ) 



k B T 



■In 



CO 



(4) 



where z is the valence and cq the bulk concentration 
of the counterions. In a dilute solution of small ions, 
this leads to cutoffs well below typical voltages for ICEO 
flows. For example, even if only steric effects are taken 
into account, with e.g. a = 3 A(including a solvation 
shell), then * c w 0.33V for c = 10~ 5 M and z = 1. To 
account for the obvious excess ions in PB theory, Stern 
long ago postulated a static compact monolayer of sol- 
vat ed ions 22]. This is also invoked in some models of 
ICEO flows, where a constant capacitance is added to 
model the Stern layer and/or a dielectric coating, which 
carries most of the voltage when the diffuse-layer capaci- 
tance ([3]) diverges. However the voltage drops applied in 
ICEO make it unrealistic that a monolayer could with- 
stand most of the drop, and furthermore a dynamical 
model is required for a condensed layer that is built and 
destroyed as the applied field alternates. 
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FIG. 1: (a) The equilibrium distribution of counterions near 
a surface in Poisson-Boltzmann (PB) and a modified (MPB) 
theory taking into account a minimum ion spacing a for 
large applied voltages ze^olksT — 20, 30, . . . , 100, also 
zety d /hsT — 7, 10 are shown for MPB to illustrate interme- 
diate stages, (b) The voltage dependence of the differential 
capacitance Cd of the diffuse layer from @ and (JjjJ, where 
v = 2a 3 co is the bulk volume fraction of ions. 



A variety of "modified Poisson-Boltzmann equations" 
(MPB) have been proposed to describe equilibrium ion 



profiles near a charged wall [23J, |24(. To capture ion 
crowding effects across a wide range of voltages, we em- 
ploy the simplest possible MPB model of Bikerman 25 1 



and others [231 . 1261 ] , which is a continuum approximation 
of the entropy of ions on a lattice of size a. As shown in 
Fig. [TJa), when a large voltage is applied, the counterion 
concentration exhibits a smooth transition from an outer 
PB profile to a condensed layer at c = c max = a~ 3 near 
the surface. 



Equation Q of dilute-solution theory predicts that Cd 
diverges with ^>d, but with this model we predict the 
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FIG. 2: (a) Sketch of one period of the most studied ACEO 
pump geometry [53], indicating slip (arrows) and induced 
diffuse-layer charge, and (b) the (dimensionless) flow rate 
versus frequency for different models. In the low-voltage 
limit V <C fcsT/e = 25 mV, all models predict a sin- 
gle peak (solid line). For a typical experimental voltage, 
V = WOksT/e — 2.5 V, Poisson-Boltzmann theory breaks 
downs and its capacitance ([3| shifts the flow to very low fre- 
quency (dot-dashed line) , but accounting for steric effects © 
with v = 0.001 (dashed line) reduces the shift and predicts 
high frequency flow reversal, similar to experiments [ill ]. 
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where v = 2a 3 cq is the bulk volume fraction of ions. As 
shown in Fig. QJb) , the capacitance reaches a maximum 
near the critical voltage ty c and then decreases at larger 
voltages because the effective capacitor width grows due 
to steric effects. 

This decrease of diffuse-layer capacitance at large volt- 
ages is robust to variations in the model and has major 
consequences for nonlinear electrokinetics. For example, 
it provides a simple explanation for the flow reversal seen 
in ACEO experiments, without invoking Faradaic reac- 
tions. As shown in Fig. [21 numerical simulations of a well 



studied planar pum p g eometry [111 . 1271 ] with the standard 
linear model 0, 0, HlJ [28j] predict a single peak in flow 
rate versus frequency at all voltages. If the nonlinear 
PB capacitance ([3]) is included then the peak is re- 
duced and shifts to much lower frequency (contrary to 
experiments), due to slower charging dynamics at large 
voltage |17| . The MPB capacitance with steric ef- 
fects (O reduces the peak shift and introduces flow re- 



versal similar to experiments, albeit with a large value 
of a = 0.0005cq , perhaps due to the underprediction of 
liquid steric effects by the lattice approximation [3]. 

The physical mechanism for flow reversal is simple: At 
low voltage, the pumping direction is set by the larger 
electrode, but at large voltages, since the more highly 
charged smaller electrode has its "RC" charging time re- 
duced by steric effects, it "wins" in pumping against the 
larger electrode at higher frequency. Perhaps a similar 
effect is responsible for the double-peaked structure at 
high voltage with non-planar stepped electrodes [HI]. 

Induced viscosity increase - The strong decay of ICEO 
flow with inceasing concentration suggests a dramatic in- 
crease in the viscosity of a highly charged diffuse layer. 
Classical continuum theory provides a general formula 
for the electro-osmotic mobility [3] , 
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as an integral over the potential difference \& entering the 
diffuse layer from the bulk. At large voltages, the effec- 
tive zeta potential £ (a measure of flow generated) should 
be smaller than as e decreases (due to alignment of 
water dipoles) and r\ increases (due to viscoelectric ef- 
fects) within the diffuse layer. 

Focusing on the viscoelectric effect in water, Lyklema 
and Overbeek [i^, [3(3] first derived a modified slip for- 
mula by assuming r\ oc E 2 in PB theory and predicted 
that b(^D) saturates at a constant value, which decays 
with increasing cq. The saturation, however, relies on 
the unphysical divergence of the counterion concentra- 
tion (and thus E) in PB theory; with MPB steric effects 
included, it can be shown that the mobility does not sat- 
urate, and the viscoelectric effect is not strong enough to 
describe ICEO experiments. 

In a very crude attempt, we adopt the MPB theory 
above and postulate that e/rj diverges as the counterion 
density (equivalent to the charge density p) approaches 
the steric limit c max as: 
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Vb 



l 



- a 3 p) 
Vb 
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This (arbitrary) choice leads to a simple formula for the 
effective zeta potential, 



C = ^d- sgn(* D ) log 



1 + 4a 3 c sinh 



ze^i 



2k B T 

( 8 ) 

whichreduces to © (C ~ *d) for < *c but satu- 
rates C ~ for ^>d 3> ^c- From Q, we recover the ex- 
perimentally observed scaling with concentration in the 
large voltage limit, u oc log(c c /co), with c c = c max = a~ 3 - 
The mobility b — EbC/Vh from is also sensitive to the 
solution chemistry, through a, z, and Co, unlike the clas- 
sical formula © valid at low voltages. The saturation 
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Choosing a "steric" value a ~ 1 — 4 A in (O shifts the 
flow to too low frequency in dilute solutions (as in PB 
theory [13]) and overestimates the concentration scale 
for ICEO flow suppression in ©. More realistic MPB 
models, which predict stronger steric effects may 
improve the fit, but correlation effects on ICEO flow may 
also need to be described. 

Of course, our models are over-simplified, but remem- 
ber that the challenge is to describe ICEO flow over 
more than three decades of diffuse-layer voltage from 
ksT/e = 25 mV to ps 10 V. The upper limit corresponds 
to a new regime for the theory of electro-osmosis, where 
countcrions are condensed near a highly charged surface. 
Nanoscale experiments and atomic-level simulations will 
be crucial to further develop the theory. 



FIG. 3: Effective zeta potential £ = u s rjb/E t eb (slip velocity 
per tangential electric field) versus diffuse-layer voltage ^ d at 
different bulk concentrations Co for the modified slip formula 
(|8|) postulating a viscoelectric effect related to crowding of 
ions in the diffuse layer at large voltages. 
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of C, also implies that the scaling of ICEO flows changes 
from quadratic, u oc E 2 or V 2 , to linear, u oc \E\ or |V| 
at large voltages. 

These predictions make the theory more realistic, but 
the experimentally observed c c = 10 mM implies a mean 
ion spacing of a = 4.4 run (roughly 40 atomic diame- 
ters) for the divergence of e/ij. How could we explain 
such a large value? Perhaps at large voltages, counteri- 
ons condense into a sort of Wigner crystal, which resists 
shear due to strong electrostatic correlations (in addi- 
tion to viscoelectric effects in the solvent). Indeed, the 
mean-field approximation breaks down when ion spac- 
ings approach the Bjerrum length, Ib — {ze) 2 / AiteksT , 
which is 7 A for bulk water and monovalent ions 
(z = 1). If e ss O.lef, (as electrochemists infer for the 
Stern layer [22]), then a w Ib is possible, so correlation 
effects on electro-osmotic flow (which to our knowledge 
have never been studied) could be very significant at large 
voltages, even in dilute bulk solutions. 

In conclusion, we have argued that (at least) two new 
phenomena arise in electro-osmosis at large induced volt- 
ages: (i) Crowding effects decrease the differential capac- 
itance (Fig. 1) which can explain high frequency flow 
reversal in ACEO pumps (Fig. 2); (ii) viscosity increase 
upon ion crowding saturates the mobility (Fig. 3), which 
implies dependence on solution chemistry and flow decay 
with increasing concentration. Although we believe these 
predictions are robust we have not managed to combine 
our simplist models (O and ([8]) into a complete theory. 
For example, choosing a ss 1 — 4 nm in J5|) to fit the 
critical concentration c c w 0.01 — 0.05 M tends to elimi- 
nate flow reversal in Fig. 2b since the reduced mobility 
of small electrode in Fig. 2a offsets its faster charging. 
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